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ABSTRACT 

We determine the underlying shapes of spiral and elliptical galaxies in the Sloan Digital 
Sky Survey Data Release 6 from the observed distribution of projected galaxy shapes, 
taking into account the effects of dust extinction and reddening. We assume that the 
underlying shapes of spirals and ellipticals are well approximated by triaxial ellipsoids. 
The elliptical galaxy data are consistent with oblate spheroids, with a correlation 
between luminosity and ellipticity: the mean values of minor to middle axis ratios 
are 0.41 ± 0.03 for Mr « -18 ellipticals and 0.76 ± 0.04 for Mr -22.5 ellipticals. 
Ellipticals show almost no dependence of axial ratio on galaxy colour, implying a 
negligible dust optical depth. 

There is a strong variation of spiral galaxy shapes with colour indicating the 
presence of dust. The intrinsic shapes of spiral galaxies in the SDSS-DR6 are consistent 
with flat disks with a mean and dispersion of thickness to diameter ratio of (21 ± 2)%, 
and a face-on ellipticity, e, of ln(e) = —2.33 ± 0.79. Not including the effects of dust in 
the model leads to disks that are systematically rounder by up to 60%. More luminous 
spiral galaxies tend to have thicker and rounder disks than lower-luminosity spirals. 
Both elliptical and spiral galaxies tend to be rounder for larger galaxies. 

The marginalised value of the edge-on r-band dust extinction Eq in spiral galaxies 
is Eq ~ 0.45 magnitudes for galaxies of median colours, increasing to i?o = 1 magni- 
tudes for g ^ r > 0.9 and Eq — 1.9 for the luminous and most compact galaxies, with 
half-light radii < 2/i^^kpc. 
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INTRODUCTION 

The quantitative study of intrinsic galaxy shapes started 
with Hubble (1930), who measured the projected axial ra- 
tios of elliptical galaxies when classifying them into what 
would later become the Hubble sequence. Using the pro- 
jected axial ratios measured from photographic plates of 
254 spiral galaxies from the Reference Catalogue of Bright 
Galaxies (de Vaucouleurs & de Vaucouleurs 1964) , Sandage, 
Freeman & Stokes (1970) concluded that the disks of spiral 
galaxies were circular, with a disk thickness (defined as the 
ratio of disk height to diameter) of 7 = 0.25. Later esti- 
mates from photographic plate surveys were performed by 
BinggeH (1980), Benacchio & Galletta (1980), and Binney 
& de Vaucouleurs (1981), who concluded that galactic disks 
were consistent with almost circular ellipses, with a mean 
ellipticity of e = 0.1. These results, based on small samples 
of galaxies, have been superseded in recent years by much 
larger studies from CCD imaging and scans of wide-field 
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photographic surveys (Fasano & Vio 1991). Lambas, Mad- 
dox & Loveday (1992) analyzed a sample of ~ 13, 000 APM 
galaxies, and found that the distribution of ellipticities was 
well fitted by a one-sided Gaussian distribution centred on 
e = with a dispersion of = 0.13 and a mean of (e) = 0.1. 
Rix & Zaritsky (1995) studied a sample of kinematically se- 
lected face-on spiral galaxies in more detail, finding a typical 
ellipticity of e = 0.045 in the galactic disk potential. 

Spatially resolved observations of internal kinematics 
can sort out the three-dimensional shape of a galaxy (Bin- 
ney 1985; Franx et al. 1991; Statler 1994ab, Statler & Fry 
1994; Bak & Statler 2000; Statler, Lambright, & Bak 2001). 
Andersen et al. (2001) and Andersen & Bershady (2003) ap- 
plied this method to 24 largely face-on spirals and found a 
mean ellipticity of (e) — 0.076, similar to that of Rix & Zarit- 
sky. However, in both cases the selection of face-on objects 
may have introduced systematic biases in the sample. Future 
work with the SAURON spectrograph (Bacon et al. 2001; 
de Zeeuw et al. 2002) will allow detailed three-dimensional 
models to be created for a much larger number of galaxies. 

Taking advantage of the large number of galaxies with 
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high-quality photometry and shape measurements in the 
Sloan Digital Sky Survey (SDSS; York et al. 2000), Ryden 
(2004) selected a sample of spiral galaxies from the SDSS 
Data Release 1 (DRl, Abazajian et al. 2003), chosen to min- 
imise systematics due to seeing. She found that the distribu- 
tion of galactic disk ellipticities can be well fit by a Gaussian 
distribution in In e with a mean of —1.85 and a standard de- 
viation of 0.89. Vincent & Ryden (2005) extended this work 
using the SDSS Data Release 3 (Abazajian et al. 2005), and 
fit the distribution of axis ratios of both ellipticals and spi- 
rals to tria^xial models. Assuming a uniform triaxiality (i.e. 
all galaxies are either prolate, triaxial or oblate), they found 
that both spiral and elliptical distributions are consistent 
with oblate spheroids. Moreover, high luminosity elliptical 
galaxies show rounder shapes than do lower-luminosity el- 
lipticals. 

Elliptical galaxies were once believed to be axisymmet- 
ric oblate spheroids, until it was discovered that their rota- 
tion velocities were insufficient to support such a geometry 
(Bertola & Capaccioli 1975). Binney (1976) suggested that 
ellipticals could be well described by a triaxial ellipsoid but 
Davies et al. (1983) found that small ellipticals are better 
fit by oblate spheroids. This variety of intrinsic shapes for 
elliptical galaxies makes it difficult to obtain their intrinsic 
shapes using only their apparent images; when this approach 
is used on large numbers of elliptical galaxies, it is often nec- 
essary to assume a triaxiality as in Vincent & Ryden (2005) , 
or to use the misalignmeiit between the internal isophotes 
of individual elliptical galaxies as suggested by Binney & 
Merrifield (1987). The study of intrinsic shapes of spheroids 
has been recently extended to bulges in spiral galaxies by 
Mcndez-Abreu et al. (2008); bulge shapes are found to be 
consistent with a mean axial ratio in the equatorial plane of 
(B/A) =0.85. 

Astronomers as early as Holmberg (1958) realised that 
the shape distribution of spiral galaxies is affected by the 
presence of dust. Optically thick dust obscuration aligned 
in the rotational plane of spirals will cause edge-on ob- 
jects to appear systematically fainter, and thus they will 
be under-represented in magnitude-limited samples, biasing 
the estimates of intrinsic galaxy shapes. The dust extinc- 
tion of galaxies is important for understanding the true lu- 
minosities of galaxies, the distribution of ISM in galaxies, 
and the relationship between optical and infrared emission 
from galaxies (for reviews, see Davies & Burstein 1995 and 
Calzetti 2001), and studies of the brightness of galaxies as a 
function of axial ratio should allow the effects of dust to be 
quantified. Valentijn (1990) studied the shapes and surface 
brightnesses of 16, 000 galaxies from digitised photographic 
plates, and interpreted the data as indicating an optically 
thick component in disk galaxies, extending well beyond the 
apparent optical extent of the galaxy. Burstein, Haynes & 
Faber (1991) and Choloniewski (1991) however, showed that 
Valentijn's results were due in part to selection effects, and 
found that the diameters of galaxies were independent of 
inclination; see Davies et al. (1993) and Valentijn (1994) 
for further discussion of these issues. Peletier & Willmer 
(1992) expanded on the effects of selection biases with incli- 
nation, and emphasised that that the dust opacity may de- 
pend on galaxy luminosity. Tully et al. (1998), for example, 
found a 1.3 mag difference in the R band between face-on 
and edge-on luminous galaxies, but found a negligible ef- 



fect for intrinsically faint galaxies. Holwerda et al. (2005a, 
2005b) used a more direct method for obtaining the opacities 
of spiral disks, consisting of measuring the rmmber of field 
galaxies seen through galactic disks using images from the 
Hubble Space Telescope WFPC2 archival data. This method 
had previously been applied to ground-based data by many 
other authors, including Zaritsky (1994), Nelson Zaritsky 
& Cutri (1998), and Keel & White (2001). Valotto & Gio- 
vanelli (2005) followed a different approach to derive the 
dust extinction in galaxies, using the inner part of the rota- 
tion curves of spiral galaxies. 

More recently, a number of groups have studied the vari- 
ation of galaxy properties with the inclination angle with 
respect to the line-of-sight, or more directly with projected 
galaxy shapes, to draw conclusions regarding dust extinction 
in spiral galaxies. Shao et al. (2007) measured dust extinc- 
tion in spiral SDSS-DR2 (Abazajian et al. 2004) galaxies 
by studying the luminosity function of galaxies with differ- 
ent inclination angles, and using the intrinsic galaxy shapes 
as inferred from the distribution of projected axis ratios. 
They interpret the decrease in characteristic LF luminosity 
(L*) with increasing inclination as an effect of dust extinc- 
tion, where the disk optical depth is roughly proportional 
to the cosine of the inclination angle. However, they did 
not take into account the influence of dust on the projected 
shapes of galaxies. Unterborn & Ryden (2008) also study 
the variation of the luminosity function with inclination us- 
ing a subsample of ~ 78, 000 galaxies from the SDSS Data 
Release 6 (DR6, Adelman- McCarthy et al. 2008), finding 
similar results for the dependence of extinction on projected 
shape. They use this to define an extinction-unbiased sample 
of spiral galaxies for which they estimate intrinsic shapes. 
Even though their results indicate that these galaxies axe 
consistent with flattened disks as was found by previous au- 
thors (e.g. Ryden 2004), the definition of the sample makes 
it difficult to compare their results with previous estimates. 
Mailer et al. (2008) study the variations of galaxy proper- 
ties with inclination and derive extinction corrections using 
the NYU-VAGC (Blanton et al. 2005), which combines data 
from SDSS and the Two-Micron All Sky Survey (Skrutskie 
et al. 2006) . The median extinction over their whole sample 
(all morphological types) is 0.3 magnitudes in the p-band. 
Finally, Driver et al. (2007) also study the dependence of 
the luminosity function with inclination, by decoposing their 
sample of galaxies in the Millenium Galaxy Catalogue (Liske 
et al., 2003, Driver et al., 2005) into bulge and disk compo- 
nents, and are able to deduce the residual face-on attenua- 
tion. 

Dust has also been found in elliptical galaxies. Ebneter, 
Davis & Djorgovski (1988) used colour maps to find evidence 
of dust in more than 30% of their sample of elliptical galax- 
ies; ~ 2.5% of the galaxies showed evidence for a dusty disk. 
However, the amount of dust in ellipticals is rather smaller 
than in spirals. For instance, Knapp et al. (1989) found that 
an elliptical galaxy contains between 1 and 10% of the dust 
content present in a spiral galaxy of similar luminosity (see 
also Leeuw et al. 2004, Krause et al. 2003, Goudfrooij 2000). 
Far-infrared observations of elliptical galaxies by Temi et al. 
(2004) also place constraints on the mass of dust in ellipti- 
cals in the range M^uat = 10^ — 10^ /i~^Mo, where h is the 
Hubble constant in units of lOOkms"^ Mpc~^. This mass 
is ~ 10~® of the stellar mass, a much smaller fraction than 
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seen in spiral galaxies, where the fraction is of order 5 x 10~ 
(see for instance Stevens, Amure & Gear, 2005). 

This paper will use the SDSS DR6 to derive the intrin- 
sic three-dimensional shapes of spiral and elliptical galaxies 
using the apparent photometric shapes of galaxies. We will 
include an in-depth analysis of the effects of dust on the 
distribution of apparent shapes of spiral galaxies, which will 
explain most of the trend seen in the distribution of spi- 
ral shapes with luminosity and colour. Our analysis weights 
galaxies by the inverse of the volume out to which they can 
be seen, thus simulating a volume-limited catalogue, and 
allowing us to use the full sample of galaxies available in 
the SDSS DR6. Furthermore, the large number of galaxies 
present in this sample allows us to study the dependence of 
the intrinsic shapes of galaxies with luminosity, colour and 
physical size. Throughout this paper we assume a standard 
ACDM cosmology, with matter density parameter Qrn = 0.3 
and a cosmological constant corresponding to flA = 0.7. 

This paper is organised as follows. In Section |2] we will 
briefly describe the SDSS DR6 galaxies, and the parame- 
ters we consider when measuring the distribution of shapes. 
Section [3] explains our methodology, including our model for 
the effects of dust on the observed distribution of axis ratios. 
Section |4] shows our results, and Section [5] summarises the 
main conclusions drawn from this work. 



2 THE SDSS GALAXY SAMPLE 

We select the ~ 585, 000 galaxies from the r < 17.77 
magnitude-limited main spectroscopic galaxy sample of the 
SDSS (Strauss et al. 2002) from the SDSS DR6 (Adelman- 
McCarthy et al. 2008). The SDSS imaging data consist 
of CCD imaging data in five photometric bands (ugriz, 
Fukugita et al. 1996), taken with a drift-scan camera (Gunn 
et al. 1998) on a dedicated wide-field 2.5m telescope (Gunn 
et al. 2006). The properties of all detected objects in the 
images are measured (Lupton et al. 1999; Stoughton et al. 
2002), and are calibrated astrometrically (Pier et al. 2003) 
and photometrically (Smith et al. 2002; Ivezic et al. 2004; 
Tucker et al. 2006). We K-correct the galaxy magnitudes us- 
ing V3.2 of the code described in Blanton & Roweis (2006). 

The image of each galaxy in the SDSS sample is fit to 
two-dimensional models of a de Vaucouleurs (1948) surface 
profile and an exponential profile, each convolved with the 
PSF of the image. This fitting procedure provides a measure- 
ment of the model axial ratios (6/a) of each galaxy image 
in a way robust to seeing, as well as the effective radius and 
position angle of the galaxy. The SDSS image pipeline also 
fits an ellipse to the 25 mag/arcsec^ isophote of each galaxy, 
and determines the so-called adaptive moments (Bernstein 
& Jarvis 2002); while the axis ratios via these statistics are 
generally in good agreement with those from the model fits, 
they are affected by seeing, and thus tend to yield system- 
atically rounder shapes than do the model fits. 

We will infer the three-dimensional shapes of spiral and 
elliptical galaxies separately. Park & Choi (2005) presented 
a very accurate way to determine SDSS galaxy morpholo- 
gies using colour gradients that could be used to separate 
the DR6 catalogue into spiral and elliptical galaxies, but 
this would require analyzing the images of each individual 
galaxy separately. We use an alternative method: in fitting 



the exponential and de Vaucouleurs models, the SDSS imag- 
ing pipeline also asks for the best linear combination of these 
models (Abazajian et al. 2004) , as quantified by the param- 
eter fracDeV. We use this parameter to distinguish spiral 
galaxies (fracDeV < 0.8) from ellipticals (fracDeV ^ 0.8). It 
is not possible using these techniques to separate out lentic- 
ular or 50 galaxies; the definition of this morphological type 
is difficult using photometric data, out to the redshifts ex- 
plored in this work. The axis ratios from the exponential 
and de Vaucouleurs models are in excellent agreement, in- 
dependent of the value of fracDeV, with a scatter of about 
0.05 around the identity line, but we adopt the axis ratios 
from the exponential fit when fracDeV < 0.8, and the de 
Vaucouleurs parameters otherwise. All our analyses are car- 
ried out using model DeVaucouleurs or exponential r-band 
magnitudes and g — r colours, depending on the galaxy type. 



3 THE INTRINSIC SHAPES OF GALAXIES 

In this section we will measure the distribution of projected 
axis ratios and present the model that will allow us to infer 
their intrinsic shapes. 

3.1 Distributions of projected axis ratios of 
elliptical and spiral galaxies 

Figure [T] shows distributions of spiral galaxies on the left 
panels, and of elliptical galaxies on the right panels. The top 
panels show distributions of galaxy redshifts, and the middle 
and lower panels the distributions of projected axis ratios. In 
addition, the thick solid lines in the top and middle panels 
show the distributions of redshifts and projected axis ratios 
for the full sample of galaxies in the SDSS-DR6. In these 
panels, the thin lines illustrate the variation of redshift and 
projected axis ratio distributions as the typical galaxy lu- 
minosity of the sample is increased. More luminous galaxies 
tend to show rounder apparent shapes (axial ratios closer 
to unity), suggesting that the intrinsic shapes of galaxies 
are a function of luminosity. This is a fiux-limited sample, 
and thus has a strong correlation between luminosity and 
redshift, as shown in the top panel. With this in mind, we 
simulate a volume-limited measurement of apparent shapes 
by simply weighting each galaxy by 1/Vmax, where Vmax is 
the volume corresponding to the maximum distance out to 
which a galaxy of a given apparent magnitude enters the 
flux-limited catalogue, taking into account K-correctioni{3. 
As a result of these procedures, our final samples of spirals 
and ellipticals contain a total of 282203 and 303390 galaxies, 
respectively. There is an additional effect for spiral galaxies, 
whereby the internal extinction will cause edge-on objects 
to appear fainter than equally luminous face-on objects. We 
thus first volume-limited ignoring this effect, determined the 
inclination dependence of the extinction as described in ij |3.2l 
below, and then redefined our volume-limited sample taking 
this extinction into account and repeated our analysis. We 

^ The advantage of using the 1/Vmax weight is the larger sample 
size and the increased range of luminosities that can be explored; 
using a true volume-limited sample would restrict our analysis to 
a sample almost ~ 10 times smaller composed only by intrinsically 
bright galaxies. 
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Figure 1. Left panels show distributions obtained for spiral galaxies, right panels for elliptical galaxies. Top panels: distribution of 
galaxy rodshifts, normalised so that the area under each curve is unity. The thick solid line corresponds to all the galaxies in the SDSS- 
DR6. The thin lines show the distribution of axis ratios for galaxies in different bins of absolute magnitude as indicated in the key, and 
for spiral and elliptical galaxies separately. Middle panels: normalised distribution of axis ratios for the same samples of galaxies as in 
the top panels; the thick line corresponds to the full sample of galaxies. Bottom panels: Distribution of axis ratios, summed over all 
luminosities. The solid lines show the results when a 1 IVmax weight is applied to each galaxy. The dotted lines show the results with no 
weighting. Errors are calculated using the Jack-knife technique. 



found that the inferred extinction differed by only 0.1 mag 
between the two analyses, so we did not iterate further. The 
solid lines in the bottom panels of Figure [1] show the re- 
sulting distributions of axial ratios for spiral and elliptical 
galaxies. Not surprisingly, the spiral galaxy distribution is 
skewed toward lower axial ratios (hja values) than are the 
ellipticals, which are rounder with hja closer to 1. 

The distribution becomes flatter when using the l/V^ax 
weighting. In particular, the weighted distribution of hja 
values for spiral galaxies is qualitatively similar to that of 
Ryden (2004), who found a fiat distribution over a wide 
range of hja values from fe/a = 0.2 to 0.7 for a volume lim- 
ited sample of SDSS-DR2 spiral galaxies. We use l/Knax 
weighting in all the analyses that follow. Galaxies fainter 
than Mr — 5 logjQ(/i) = — 17 have very small values of Vma^, 
and thus tend to dominate the noise of the estimate of the 
distribution function; we therefore drop such low-luminosity 
objects in what follows. Throughout this paper, error-bars 
were calculated using the Jack-knife method. 

We model galaxies as triaxial ellipsoids of major axis 
A, middle axis B, and minor axis C, parameterised by two 
axis ratios, C/B and B jA, and we will determine the distri- 



bution of axis ratios of the spiral and elliptical populations 
separately. Following Ryden (2004), we assume that the dis- 
tribution of 1 — C/-B of the three-dimensional structure can 
be approximated by a Gaussian with mean 7 (this parameter 
is related to ^.y in Ryden 2004, via 7 = 1 — /x-y) and standard 
deviation a~f. We also assume that there is a log-normal dis- 
tribution in the quantity e = log(l — -B/A), with mean [i 
and dispersion cr. Larger values of 7 and [s, correspond to 
more elliptical objects in the B — C and A — B planes, re- 
spectively. Given values of axis ratios drawn from these dis- 
tributions, and a random viewing angle {9,(f)), we compute 
the resulting apparent axis ratio, b/a (Binney 1985), using 
equations 12 — 15 from Ryden (2004). Repeating this mul- 
tiple times gives a model distribution (A^modci (&/«)) which 
can be compared directly to our measured volume-weighted 
distributions {N{b/a)). 

This model assumes that the dust extinction in galax- 
ies is independent of the viewing angle. However, we can 
test for the effects of dust by exploring the dependence of 
the shape distribution on absolute magnitude and colour. 
The top-left panel of Figure [2] shows the median h/a of 
spiral and elliptical galaxies as a function of r-band abso- 
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Figure 2. Top-left: median axis ratios of spiral (short dashed) and elliptical galaxies (long dashed) as a function of r-band absolute 
magnitude; errors are calculated using the Jack-knife technique. The middle-left panel shows the distributions of absolute magnitudes 
(luminosity functions) as short and long dashed lines. The filled circles show the best Schechter fit to SDSS galaxies from Blanton et al. 
(2003a) normalised so that the area under the curve is equal to one. The bottom-left panel shows the luminosity functions for all spiral 
galaxies, and for the subset of face-on spiral galaxies (b/a > 0.8) which we use as our unextincted luminosity function. Right: Median 
axis ratio and absolute magnitude as a function of g — r colour for spirals and ellipticals separately. We do not show results for the full 
sample to improve clarity. 



lute magnitude. Both spiral and elliptical galaxies tend to 
be rounder at higher luminosities. Spiral galaxies show a 
change of A(fe/a) ~ 0.2 between absolute magnitude values 
of Mr — 51ogio(/'-) ~ —18.5 and —22.5. The variation in b/a 
for ellipticals is smaller, A{b/a) ~ 0.1. 

The top-right panel shows the median b/a as a function 
oi g — r colour. Red spiral galaxies show systematically larger 
axis ratios than do blue spirals. We interpret the colour and 
luminosity dependence of spirals as due to the reddening and 
dimming effects of dust, which becomes more prominent for 
edge-on systems. Elliptical galaxies show only a mild vari- 
ation in shape with absolute magnitude, and no significant 
effect in colour. Thus, not unexpectedly, we see no evidence 
for a dust layer aligned with the principal plane of elliptical 
galaxies, and assume that the trend in the apparent shape 
of ellipticals with luminosity corresponds to a real variation 
of their intrinsic shapes. We study the dependence of pro- 
jected and intrinsic properties of spiral and elliptical galaxies 
in more detail in Section |4l 



The middle panels of figure [2] show the luminosity 
and colour distribution functions (left and right panels, 
respectively), both estimated using the 1/Vmax estimator, 
for spiral and elliptical galaxies. The filled circles in the 
middle-left panel give the r-band luminosity function esti- 
mate from Blanton et al. (2003a), with Schechter parameters 
Af* -51ogio(/i) = -20.44 and a = —1.05; our results for the 
full sample (black solid lines) are in excellent agreement. As 
the efi'ects of dust are more severe for edge-on objects, we 
also calculate the luminosity and colour functions for face- 
on objects {b/a > 0.8) only. These are shown as solid lines 
for the spiral galaxies in the bottom panels of this figure; 
the dashed lines show the luminosity and colour functions 
for all the spiral galaxies for comparison. The break in the 
face-on luminosity function is shifted towards higher lumi- 
nosities, as expected for the sample of galaxies least affected 
by dust. This measurement of the "unextincted" luminos- 
ity function is in agreement with results from Shao et al. 
(2007), Unterborn & Ryden (2008) and, taking into account 



6 Nelson D. Padilla and Michael A. Strauss 



J 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [ 1 1 1 L 


J 










C) 


C) 






J_L 




. , 1 . 


, . 1 , . , , 1 1 





-18 -19 -20 -21 -22 

J I I I I I I I I I I I I I I I I I I I I I I I I I I 



0.5 - 




- 



Figure 3. Dependence of the best-fitting parameters on galeixy 
luminosity, for elliptical galaxies. Squares represent the best 
fit parameters from the marginalised, one-parameter probabili- 
ties. Circles correspond to the best fit parameters in the four- 

dirncnsional parameter space. Dotted lines and shaded areas in- 
dicate the best fit parameters for the full sample of elliptical galax- 
ies. Top panel: variations in the typical fj, = median log(l — B/A) 
axis ratio; the error bars and shaded areas do not correspond to 
uncertainties in the parameters but indicate the best fitting width 
a for the Gaussian distribution of // values used in the model. 
Bottom panel: same as top panel for 7 = median(l — C/B). 
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Figure 4. Comparison between the best fit model axis ratio 
distributions and the actual measured distributions from elliptical 
galaxies in the SDSS DR6. Errors are calculated using the Jack- 
knife technique. 



the range allowed by the analysis, with Mailer et al. (2008). 
On the other hand, the colour function of spiral galaxies is 
shifted to the red due to the effects of dust. These estimates 
of unextincted luminosity and unreddened colour functions 
will be needed when we model the effects of dust, a subject 
to which we now turn. 

3.2 Modelling the effects of dust 

In this section, we develop a simple model for the impact of 
a planar distribution of dust on the distribution of appar- 
ent axis ratios of spiral galaxies. We assume no correlation 
between the dust column and the physical diameter of the 
galaxy, an assumption we will justify a posteriori. While Un- 
tcrborn & Ryden (2008) use the inclination dependence of 
the luminosity function to define a sample of galaxies not 
affected by dust (i.e. not biased towards face-on objects), 
our shape fitting solves for the dust effects solf-consistcntly. 

We follow the following steps to produce the predicted 
distribution of projected axis ratios, given our assumed dis- 
tribution of axis ratios and our measured luminosity and 
colour distributions: 

(i) We assume that the amount of extinction and red- 
dening are roughly proportional to the path length of the 
light through the galaxy. Therefore, we expect a minimum 
extinction when a galaxy is seen face-on, and an increas- 
ing extinction as the line-of-sight approaches the plane of 
the galactic disk. Similarly, Shao et al. (2007), Unterborn & 
Ryden (2008) and Mailer et al. (2008) all find that the op- 
tical depth increases monotonically with inclination angle. 
The following parametrisation of the angle dependence of 
dust extinction is not intended as a physical model, but as 
a heuristic guess for the scaling. Consider an oblate trieixial 
galaxy with axis ratios given by x = BjA and y — C/B. 
The total dust extinction as a function of inclination 6 in 
our model is 



Eo{l + y — cos 6*), if cos ^ > y 
Eo , if cos 9 < y 



(1) 



m = 



(2) 



where Eq is the edge-on extinction in magnitudes in a given 
band and y is the galaxy height to diameter ratio extracted 
from a distribution of mean 7 and width a-^ . The same can 
be assumed for the dust reddening, 

Jio(l + y — cosO), if cosO > y 
Jio,if COS0 < y 

where Ro is the edge-on reddening in magnitudes. In the 

optically thin case, we can tie Ro to the extinction via Eq = 
2.77R0, as is appropriate for the r-band and g — r colour. 
We find that our model results are not strongly dependent 
on this parameter and therefore can be applied to either the 

optically thin or thick case. 

(ii) We produce an extincted LF defined by 

4>E{M,e)=4>{M + E{e)) 

where (f){M) is the unextincted luminosity function calcu- 
lated using only face on galax;ies. 

We define the ratio, fs, between the number of ob- 
served (extincted) and intrinsic galaxies of a given luminos- 
ity /is(M) = 4>E{M)/(j>(M). Similarly, we define /fl(ff-r) as 
the ratio between the underlying and reddened distributions 
of galaxy colours. 
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Figure 5. Comparison between the best fit model axis ratio distributions and the actual measured distributions from elliptical galaxies 
in the SDSS DR6 in bins of absolute magnitude. Errors are calculated using the Jack-knife technique. 



(iii) We calculate the ratio of the number of galaxies seen 
at inchnation 6 to the number expected without extinction, 
by multiplying the effects of reddening and extinction to- 
gether, 



rjrjE{M)fn{C)MM)MC)W{C,M)dCdM 
/-I MM)MC)W{C, M)dMdC 



(3) 



where C — g — r, and the function W contains the correla- 
tion between colour and Mr- We assume that W is Gaussian 
with mean and dispersion extracted directly from the data; 
this correlation is compatible with the results shown in Fig- 
ures 11 and 12 of Blanton et al. (2003b). The sub-index 
s indicates that the luminosity and colour functions corre- 
spond to a particular subsample of galaxies; these subsam- 
ples are defined using sharp cuts in the allowed luminosity 
and colour ranges. This indicates that depends not only 
on the amount of extinction and reddening, but also on the 
range of luminosities and colours present in each subsample 
of galaxies. Note that Eq. [3] ignores the presence of large- 
scale structure, which is justified given the large solid angle 
of the SDSS sample. 

(iv) We then construct the model distribution of apparent 
axis ratios A'inodci(&/i) as described in Section [S] Rather 
than selecting the cosine of the viewing angle from a flat 



distribution, we select it from ^p{9) as given in equation (|3]). 
In the case where Eo = and Ro — 0, tp{0) — 1. 

In general, the effect of dust extinction and reddening is 
to decrease the number of galaxies seen edge-on relative to 
those that are face-on. This decrease depends strongly on the 
luminosity and colour functions, as well as on the selected 
range of luminosities and colours. The above process gives 
a prediction for the observed projected axis ratio distribu- 
tion for a given set of triaxial galaxy parameter distribution 
functions and dust properties. Note that for a given viewing 
angle, the model states that the dust affects the likelihood 
that a galaxy would enter the sample at that viewing angle, 
but does not affect the observed projected ax;is ratio. This 
is a good approximation as long as the dust is smoothly 
distributed within individual galaxies. 

In the following section, we will constrain these param- 
eters by fitting these predictions to the observed axis ratio 
distribution. 



3.3 Parameter fitting 

We aim to constrain the parameters /i, ct, 7 and which de- 
scribe the intrinsic shapes of galaxies by fitting the observed 
axis ratio distribution. Spirals and ellipticals have intrinsi- 
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Table 1. Initial parameter grid for Ellipticals 



Parameter 


Min. Value 


Max. Value 


Number of steps 




-4.05 


-0.05 


21 


IT 


0.4 


3 


18 


T 


0.09 


0.99 


19 




0.01 


0.36 


12 



cally different shapes, and we fit to the two separately. For 
spirals, we also include the effects of dust via the extinction 
parameter Eo which in turn defines the reddening Ro. 

We define a grid in parameter space p (four parameters 
for ellipticals, five for spirals) . Using the parameters of each 
grid point pi, we generate random three-dimensional axis 
ratios from the assumed distribution, observed at a random 
orientation (modulated by the effects of dust, Eis described 
above). We then generate ten independent model distribu- 
tions of projected axis ratios, each containing as many galax- 
ies as the sample of galaxies under analysis. We take the 
average of these ten distributions as the final model dis- 
tribution, Afmodei(&/a, {p}i)i S'lid use the jack-knife errors, 
crjack-knife(6/a), obtained from the observed distribution to 
define a between the real data, N{b/a), and the model, 

' \ 2 

2/r 1 N I N„,odci{b/a,{p}i) ~ N{b/a) 

X {{P\i) = 2^ 

b/a bins 



f7'jack-knifc(V Cl) 



(4) 



The best fit parameters correspond to the minimum value of 
throughout the parameter grid. Throughout this analysis, 
we use a bin size of A{b/a) —0.1 in presenting the observed 
and model distributions ol b/a and when calculating x^- 



4 RESULTS 

4.1 Elliptical galaxies 

The grid of parameters we used for elliptical galaxies is 
shown in Table [1] The grids include the full range of val- 
ues suggested for these parameters in the literature. The 
number of steps for each parameter, shown in the fourth 
column, corresponds to a coarse initial grid; once the best 
fit parameters are found we re-do the analysis on a finer grid 
with twice the resolution centred on the best fit parameters. 
We repeat this refinement twice. Given the lack of colour 
dependence of the axis ratio distribution for ellipticals, we 
have assumed no dust. In order to explore the dependence 
of galaxy shape on luminosity, we divide the ellipticals by 
r-band absolute magnitude in four bins, with boundaries 
given by Mr - 51ogio(/i) = -24, -21, -20, -19 and -17. 

We calculate the marginalised one-parameter likeli- 
hoods (normalised to the maximum likelihood) resulting 
from fitting the observed N{b/a). The resulting best-fit pa- 
rameters are shown in Figure |3] The "error bars" correspond 
to the best fit widths, o and a-^, of the distributions of 
log(l — B/A) and 1 — C/B. As we guessed from the pro- 
jected axis ratios, more luminous elliptical galaxies are con- 
sistent with a rounder underlying shape: the mean axis ratio 
7 changes from 0.6 to 0.2 with increasing luminosity (see Ta- 
ble |2] for values of 7 and estimated errors not shown in the 
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. 99% confidence 
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Figure 6. Marginalised likelihood contours in the ^i-Eq plane, 
for the full sample of spiral galaxies, corresponding to 1 — , 2— and 
3 — (T confidence levels (shown as solid, dashed and dotted lines, 
respectively). 




Figure 7. Comparison between the best fit model axis ratio 
distributions (solid line) and the actual measured distributions 
from the full sample of spiral galaxies in the SDSS DR6 (open 
symbols with error bars). Errors are calculated using the Jack- 
knife method. 



figure). The quantity ^ is much more constant, with the im- 
plied B/A three-dimensional axis ratio varying only slightly, 
from ~ 0.95 to ~ 0.92. 

In this figure, squares show the best-fit parameters cor- 
responding to the marginalised one-parameter maximum 
likelihoods. The open symbols show the best-fit parameters 
as obtained from the full parameter space. In most cases, 
the two estimates agree reasonably well. Table [2] shows the 
values of the best fit parameters for the four subsamples of 
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elliptical galaxies, as well as the maximum likelihood value. 
The reduced is ^ 2, indicative of a good fit, except for 
the highest luminosity subsample and the total sample. In 
these case, our simple model does not give a statistically rig- 
orous good fit, but it still follows the shape of the observed 
distribution quite well. Figure |3] directly compares the ob- 
served axis ratio distribution with the model fits for the full 
sample of elliptical galaxies. Figure [5] shows the comparison 
in each bin of absolute magnitude; the agreement between 
the model and observed distributions, while not perfect in 
every case, is impressive, and thus the four-parameter model 
is adequate to describe the real distributions of axis ratios 
for the galaxy luminosities explored here. 

Vincent & Ryden (2005) used slightly different analysis 
techniques and a sample a quarter of the size the one used 
here, but they found results in very good agreement with our 
own, namely that low luminosity elliptical galaxies are more 
consistent with prolate spheroids than are high-luminosity 
ellipticals. 

The analysis of the marginalised two-parameter likeli- 
hoods indicates that there are little or no degeneracies be- 
tween the parameters used to fit the projected shapes of 
elliptical galaxies. 

4.2 Spiral galaxies 

The distribution of spiral galaxy axis ratios depends both 
on absolute magnitude and colour (Figure [2}, which we in- 
terpret as the effect of dust in the rotational plane of the 
galaxies. We model this as described in § 13.21 giving a five- 
parameter model. We use the grid of parameters from Table 
[3] to find the model quantities that best reproduce the dis- 
tribution of projected spiral galaxy shapes. 

The location of the minimum is given in Table U] 
covariance between the parameters in our model explains 
the discrepancies with the marginalised values. As can be 
seen, in most cases the values are small, which indicates 
an excellent agreement between the model and the data. 

This corresponds to a likelihood more than two orders 
of magnitude higher than the best fit found by Ryden (2004) 
for the spiral galaxies in the Data Release 1 of the SDSS. 
Our inclusion of a dust model is partly responsible for the 
improvement in the agreement between model and data; the 
best fit model with no dust to the full SDSS-DR6 sample of 
spiral galaxies is characterised by x^/dof = 1.14, somewhat 
higher than for the model with dust, for which x^/dof = 
0.41. The sample used by Ryden (2004) is also a factor of 
~ 25 smaller than our full sample of spiral galaxies, and is 
therefore prone to higher sample variance. 

It should be noted that dust extinction and the param- 
eter n are somewhat degenerate (none of the other parame- 
ter pairs show appreciable degeneracy). Figure |6] shows the 
marginalised /i vs. Eq likelihood contours for the full sam- 
ple of spiral galaxies. Interestingly, the axis ratio distribu- 
tion alone allows the detection of extinction at only slightly 
better than 1 a. However, the strong relationship between 
colour and axis ratios shown in figure [J] and the results by 
Shao et al. (2007), Unterborn & Ryden (2008) and Mailer et 
al. (2008), strongly indicate the presence of dust extinction 
in spirals. 

Figure [T] compares the observed spiral axis ratio distri- 
bution with that from our model; the two are in excellent 



Table 3. Initial parameter grid for Spirals 



Parameter 


Min. Value 


Max. Value 


Number of steps 


Eo 


0.0 


2.7 


14 




-4.05 


-0.05 


21 


cr 


0.4 


3.0 


18 


7 


0.09 


0.99 


19 




0.01 


0.36 


12 



agreement. This figure uses the full spiral sample; we have 
assumed that the dust and shape properties of spirals are 
independent of luminosity. We test that assumption in Fig- 
ure [S] which shows the observed b/a distribution for spirals 
in different luminosity ranges, as well as the model predic- 
tion (dashed lines) . The luminosity dependence in this model 
comes about solely from the angular selection function from 
equation ([3]) through the different ranges of luminosity that 
define each sample. There is excellent agreement between 
the low luminosity spirals and the model; indeed, the low- 
luminosity subsample has a b/a distribution very close to 
that of the full sample, which is a consequence of our 1/Vmax 
weighting. This agreement progressively degrades as we go 
towards higher luminosity. Indeed, as we saw for ellipticals, 
the high- luminosity spirals tend to have larger axis ratios 
(i.e., to be more round) than do low-luminosity objects. 

There are two ways we might model this effect. A higher 
dust column for more luminous galaxies could preferentially 
remove edge-on objects from the high-luminosity bin, as sug- 
gested by Huizinga & van Albada (1992). Alternatively, as in 
ellipticals, there could be a direct correlation between three- 
dimensional shape and luminosity for spirals (e.g., Giovanelli 
et al. 1995). Indeed, there is a strong correlation between 
Hubble type and luminosity whereby high-luminosity spi- 
rals tend to be early types with large bulges (Roberts & 
Haynes 1994, Tasca & White 2005). This can explain the 
low number of small b/a objects at such luminosities. 

We test these two options by fitting our model sepa- 
rately to galaxies in each range of absolute magnitude shown 
in Figure |8] the solid lines show the best fit distributions 
with Eo as a free parameter, and the dashed lines the best 
fits without dust {Eo = 0). Figure |5] shows the dependence 
of the best-fit parameters as a function of absolute mag- 
nitude. The marginalised estimates are in good agreement 
and are not shown. Filled symbols show results with no dust, 
and open symbols show the results when dust extinction is 
taken into account. As in Figure O error bars indicate the 
1 — cr width of the distributions in jj. and 7. The sub-boxes 
in the middle and bottom panels show the model distribu- 
tion functions of x = B/A and y = C/A for each subsample 
in order of increasing sample luminosity from left to right; 
the model with dust is shown as solid lines, without dust as 
dashed lines. As can be seen, the main change in the inferred 
parameters is the value of which is systematically lower 
when dust extinction is not considered (this is a 1 — ct effect 
for the two most luminous subsamples). This can also be 
seen to a lesser extent for the parameter 7. We find that the 
roundness of disks does not change significantly with galaxy 
luminosity. The more luminous spiral galaxies show higher 
height to diameter ratios probably due to the presence of 
larger galactic bulges. 
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Table 2. Best fit model parameters for elliptical galaxies. 



Par. 


Sample 1 

-17 > Mr > -19 


Sample 2 

-19 > Mj. > -20 


Sample 3 

-20 > > —21 


Sample 4 

-21 > Mr > -24 


All 


M 


-2.85 ±0.30 


-3.05 ± 0.20 


-2.75 ±0.10 


-3.85 ±0.15 


-2.2 ±0.1 


cr 


1.15 ±0.35 


1.00 ±0.05 


2.60 ±0.15 


2.35 ±0.20 


1.4 ±0.1 


7 


0.41 ±0.03 


0.36 ±0.06 


0.56 ±0.02 


0.76 ± 0.04 


0.57 ±0.06 




0.17 ±0.03 


0.21 ±0.02 


0.25 ±0.02 


0.17 ±0.01 


0.21 ±0.02 


xVdof 


0.41 


2.0 


1.72 


7.2 


6.8 
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Figure 8. Comparison between the best fit model axis ratio distribution to the full sample of spiral galaxies (dashed line), and the 
actual measured distributions from samples of spiral galaxies with different luminosities in the SDSS DR6 (open symbols with error 
bars). Errors are calculated using the Jack-knife method. Solid lines show the results of fitting the model to each individual subsample 
including dust extinction. The dotted lines show the best fits when not including dust. 



The inferred intrinsic shapes of spiral galaxies are ro- 
bust to the effects of dust. The amount of extinction Eo ap- 
pears to be independent of luminosity, but the uncertainty 
in this parameter is large. The lower values of 7 ~ 0.6 char- 
acterising bright galaxies results in high values of projected 
b/a regardless of the viewing angle. This, in conjunction with 
the effect of extinction to reduce the number of objects seen 
edge-on, means that Eq is not very well constrained from 
this analysis. 



4.3 Axis Ratio Dependencies on Colour and Size 

We now explore the different intrinsic shapes of galaxies ac- 
cording to their g ~ r colour. This analysis may provide 
better constraints on the dust extinction in spiral galax- 
ies. Figure [10] shows the axis ratio distribution in bins of 
g — r. The dotted lines in this figure show the model using a 
fixed galaxy shape and dust corresponding to the best fit to 
the fuU spiral galaxy sample. The agreement is reasonable 
for all except the reddest subsample. Therefore, we allow 
the parameters to vary (solid lines), and the fit to the dif- 
ferent colour subsamples is improved significantly. Table [5] 
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Figure 9. Dependence of the best-fitting parameters on galaxy luminosity, for spiral galaxies. Circles correspond to the best fit 
parameters in the five-dimensional parameter space. Open symbols show results when including dust, filled circles, results with no dust. 
The horizontal lines indicate the best fitting parameters obtained from the full sample of spiral galaxies (dotted lines show the best fit 
when including dust, thick dashed with no dust). Shaded areas and thin dashed lines show the variance in the parameters that best fit 
the observed distributions of projected axis ratios. Top panel; variations in the best-fitting values of extinction (filled symbols, dashed 
line). Error-bars show the ranges of extinction values such that the likelihood is above AC = 0.36. Middle panel: variations in the typical 
/i = log(l — B/A) axis ratio; the error bars indicate the best fitting width a-y for the Gaussian distribution of 7 values used in the model. 
The corresponding distributions for each subsample are shown on the right in individual sub-panels. Bottom panel: same as middle panel 
for 7 = 1 — C/B (thick dashed lines are not shown to improve clarity) . 



Table 4. Best fit model parameters for spiral galaxies: Pull sample and dependence on 
luminosity. 



Par. 


Sample 1 

-17 > Mr > -19 


Sample 2 

-19 > Mr > -20 


Sample 3 

— 20 > Mr > -21 


Sample 4 

— 21 > Mr > —24 


All 


Eo 


0.20 ± 0.61 


0.72 ± 0.68 


0.8 ±0.7 


0.72 ± 0.49 


0.44 ± 0.24 


1^ 


-2.13 ±0.38 


-2.41 ± 0.48 


-2.17 ±0.41 


-2.17 ±0.34 


-2.33 ±0.13 


a 


0.73 ±0.19 


0.76 ±0.25 


0.70 ± 0.35 


0.79 ± 0.31 


0.79 ±0.16 


7 


0.79 ±0.02 


0.79 ± 0.03 


0.74 ±0.03 


0.62 ±0.04 


0.79 ± 0.02 




0.048 ± 0.007 


omi ± 0.007 


0.06 ±0.01 


0.11 ± 0.02 


0.0.50 ± 0.015 




0.21 


0,1!) 


l,:5l 


0. i:-! 


0. 11 
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Figure 10. Comparison between the best fit model axis ratio distribution to tlie full sample of spiral galaxies for different g — r colour 
cuts (dotted lines), and the actual measured distributions from samples of spiral galaxies with different g — r colours in the SDSS DR6 
(open symbols with error bars). Errors are calculated using the Jack-knife method. Solid lines show the results of fitting the model to 
each individual subsample including dust. 



Table 5. Best fit model parameters for spiral galaxies: dependence on galaxy 
colour. 



Par. 


Sample 1 

-0.2 < g - r < 0.6 


Sample 2 

O.G < g - r < 0.75 


Sample 3 

0.75 < g - r < 0.85 


Sample 4 

0.85 < a - I- < 1.1 


Eo 


0.4 ±0.3 


0.6 ±0.6 


0.7 ±0.5 


1.0 ±0.5 


M 


-2.13 ±0.41 


-2.77 ±0.45 


-2.45 ± 0.38 


-2.41 ±0.33 


cr 


0.7 ±0.4 


0.61 ±0.31 


0.91 ±0.36 


0.73 ±0.32 


7 


0.80 ±0.03 


0.80 ±0.03 


0.80 ±0.02 


0.79 ±0.02 




0.054 ±0.012 


0.054 ± 0.008 


0.052 ±0.006 


0.050 ± 0.005 


xVdof 


0.33 


0.17 


0.15 


0.20 



shows the best fit parameters for each colour subsample. 
Figure [TT] shows the resulting best fit parameters with and 
without dust (open and filled symbols, respectively). The 
intrinsic shapes of galaxies do not show a clear dependence 
with colour, but the inferred amount of dust extinction is 
consistent with an increase from Eo — 0.4 for blue galax- 
ies, to i?o = 1 for the reddest sample. There are only small 



differences between the recovered values of /j, and 7 when 
including dust. 

Up to this point we had assumed that there is no depen- 
dence of shape and dust extinction with the physical size of 
the galaxy. The best-fit parameter values are those appro- 
priate for galaxies of median size within each subsample. 
However, the intrinsic shape and the amount of dust ex- 
tinction of galaxies might depend on their sizes; galaxies of 
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Figure 11. Dependence of the best-fitting parameters on galaxy colour, for spiral galaxies. Panels, lines and symbols are as in Figure|9] 



diflterent sizes but similar luminosities might have difTerent 
dynamical histories, and the dust column might reasonably 
be larger in larger galaxies. 

The top panel of Figure [T^ shows the dependence of 
the median b/a on projected galaxy size as given by the 
photometric model scale-length Vgai (exponential or de Vau- 
couleurs depending on the value of the fracDeV parameter) . 
The median b/a decreases significantly for larger galaxies 
for both spirals and ellipticals. This effect is present at all 
galaxy luminosities. 

With this dependence in mind, Figure [131 shows the re- 
sults of fitting our model to samples of galaxies in bins of 
physical size at constant luminosity, where the left panels 
show the results for elliptical galaxies (parameters ^ and 
7), and the right panels for spiral galaxies (parameters /i, 
7, and Eq). The best-fit parameters are presented in Ta- 
ble [6] for faint, spiral and elliptical galaxies, and in Table 
\7\ for bright galaxies. At a given luminosity, small ellipti- 
cal galaxies have similar model parameters as those of the 
full sample of ellipticals, although at high luminosity, they 



tend to show slightly more elongated shapes (as reflected 
in the jj, parameter). Large elliptical galaxies, on the other 
hand, are more elongated than the full sample. In partic- 
ular, the Tgai > 7 kpc sample can be identified with 
elongated prolate shapes. Quantitatively, at low luminosi- 
ties {Mr = —18), the median axis ratios are B/A = 0.94 
and C/B = 0.42 for small galaxies, and B/A = 0.85 and 
C/B = 0.1 for large galaxies. We examined the SDSS im- 
ages and colours of ellipticals of the most extreme axis ratio 
{b/a < 0.1); all appeared to be correctly classified. 

Low-luminosity spiral galaxies, on the other hand, go 
from median axis ratios of B/A = 0.9 and C/B = 0.3 for 
small sizes to B/A = 0.36 and C/B = 0.12 for large galaxies. 
Thus, the disk thickness decreases with increasing size. Shao 
et al. (2007) assumed that n and a axe independent of size, 
whereas we found that most of the changes in the shape 
of galaxies with size are actually absorbed by variations in 
these parameters. However, the conclusions of Shao et al. 
are robust to this detail and indicate, as our results do, that 
larger spiral galaxies tend to have flatter disks. 
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Figure 12. Variation of projected axis ratio b/a with the physical size of galaxies. Dashed lines correspond to spiral galaxies and solid 
lines to elliptical galaxies. Each panel corresponds to a different luminosity bin. 



Table 6. Best fit model parameters: dependence on galaxy size for galaxies 
with Mr > —19. Small galaxies satisfy Vgai < 2h~^kpc, medium galaxies, 
2h~^kpc< Vgai < 7h~^kpc, and large galaxies, rgai > 7h~^kpc. 



Par. 
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Spirals 


Spirals 


Ellipticals 


Ellipticals 


Ellipticals 




Small 


Medium 


Large 


Small 


Medium 


Large 


Eo 


0.6 


0.16 


0.16 










-2.29 


-1.73 


-0.45 


-3.18 


-2.03 


-1.18 


a 


0.76 


0.64 


1.54 


0.75 


1.6 


1.6 


7 


0.71 


0.83 


0.89 


0.445 


0.325 


0.19 


(T-j, 


0.05 


0.02 


0.01 


0.17 


0.17 


0.05 



5 DISCUSSION AND CONCLUSIONS 

In this paper we have addressed the problem of reproducing 
the observed distributions of projected axis ratios of galaxies 
from the SDSS. We have introduced a number of improve- 
ments over previous works including, i) the use of larger 
samples of galaxies made possible by the introduction of an 
iterative 1/Vmax weighting scheme, ii) the inclusion of the 
effects of dust extinction on the distribution of apparent b/a 
axis ratios for spiral galaxies, iii) the analysis of dependence 



of galaxy shapes on galaxy luminosity, colour, and physical 
size. 

We developed a simple model for the effects of dust 
extinction on the distribution of apparent axis ratios and 
used it to constrain the intrinsic shapes of spiral galaxies. 
We characterise a given galaxy as a triaxial ellipsoid of axes 
A, B, and C from major to minor. The full sample of spiral 
gala^xies is characterised by the mean and standard devia- 
tion of 1 - C/B of /u = -2.33 ± 0.13 and a = 0.79 ± 0.16; 
the distribution of log(l — B/A) is modeled as a lognor- 
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Figure 13. Best-fitting parameters as a function of galaxy luminosity for diff'erent galaxy sizes (in diff'erent symbols as shown in the 
figure key). Left panels correspond to elliptical galaxies, and right panels to spiral galaxies. The top panel shows the variation of the 
extinction, Eg, middle panels show the mean value of fi, and bottom panels show the variation of the 7 parameter. 



Table 7. Best fit model parameters: dependence on galaxy size for galaxies 
with Mr < —21. Small, medium and large galaxies are selected as in table l6l 
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1.90 


0.72 


0.59 










-3.17 
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-2.13 


-2.52 


-1.12 


-3.37 


fT 


0.91 


0.58 
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0.31 


0.48 


0.66 


0.795 
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0.04 


0.14 


0.09 


0.22 
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0.17 



mal with mean 7 — 0.79 ± 0.02 and standard deviation 
cr^Y = 0.050 ± 0.015. These values are in good agreement 
with the thickness of disks derived by Ryden (2004), who 
finds 7 = 0.216, although her face-on disks are more elliptic, 
/X = —1.85, and the distribution widths are somewhat dif- 
ferent from ours. More recently, Unterborn & Ryden (2008) 
analyse the shapes of a sample of spiral galaxies which has 
been corrected for the biases introduced by dust on the dis- 
tribution of projected shapes. The resulting spiral shapes 
are characterised hy fj, — —2.56, a = 0.91, 7 — 0.216 and 



a-y — 0.067, which are roughly consistent with the results 
we present here; however, their sample selection is explic- 
itly dependent on their model for the dependence of dust on 
inclination angle, which makes it difficult to make a more 
quantitative comparison. 

We also studied variations in the intrinsic shapes of 
galaxies with galaxy luminosity (Table |4]), colour (Table [5]), 
and physical size (Tables [S] and [7|. As luminosity and colour 
are correlated, and in particular more luminous spiral galax- 
ies show larger bulges, we find rounder spiral galaxies at 
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larger luminosities and g ~r colours. At a given luminosity, 
larger galaxies tend to be flatter. 

Our quoted parameters are obtained by marginalising 
over the extinction parameter, E^. The extinction is not 
well-constrained by the axis-ratio distribution alone, but 
for our full sample, we find a value Eq = 0.44 ± 0.24 af- 
ter marginalising over other parameters, in good agreement 
with results from the literature. 

Using our results from Table |4j we determine the re- 
lation between the observed projected axis ratios of spiral 
galaxies and the inclination angles of their disks for four 
different luminosity ranges. In order to do this, we simply 
calculate {b/a)max at which the distribution of projected 
axis ratios peaks, as well as the 10% and 90% percentiles, 
for a given narrow range of viewing angles. We use the per- 
centiles to infer the range of angles that correspond to a 
given b/a ratio that can be measured from a spiral galaxy. 
Figure [14] and Table [8] show the relation between projected 
axis ratio b/a and polar viewing angle 9 for spiral galaxies 
of different luminosities. The shaded area shows the ranges 
of polar viewing angles that correspond to a given value 
of b/a for the fainter sample. As brighter galaxies tend to 
have thicker disks, the value of b/a corresponding to edge-on 
bright galaxies is higher than for fainter galaxies. Note that 
these results take into account the effects of dust extinction. 

The axis ratio distribution of elliptical galaxies shows 
no dependence on colour, suggesting that dust extinction 
is not important for this sample, and we do not include it 
in our modelling. The full sample of elliptical galaxies are 
characterised by parameters fj, = —2.2 ±0.2, a — 1.4 ±0.10, 
7 — 0.57 ± 0.06 and = 0.21 ± 0.02, which correspond 
to slightly oblate spheroids in agreement with results by 
Vincent & Ryden (2005). More luminous ellipticals tend to 
be rounder, although ellipticals are oblate at all luminosities. 
Although Vincent & Ryden (2005) found that the lowest 
luminosity ellipticals are best fitted by prolate spheroids, 
their inferred C/A axis ratio is consistent with our estimates. 

This work presented the most detailed statistical study 
of shapes of galaxies separated in spirals and ellipticals to 
date, and includes the effects of dust on the distribution of 
projected shapes of spirals in a self-consistent way. It is re- 
markable that dust can be inferred from the distribution of 
projected axis ratios alone, although only at a low statisti- 
cal significance. The dependence of intrinsic galaxy shapes 
and extinction with luminosity, colour and size expands the 
range of tests that galaxy formation models need to satisfy, 
to continue improving the modelling of the processes that 
drive the evolution of galaxies. 
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